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ABSTRACT 

The boundaries of cold dark matter halos are commonly defined to enclose a density contrast A relative to a 
reference (mean or critical) density. We argue that a more physical halo boundary choice is the radius at which 
accreted matter reaches its first orbital apocenter after turnaround. This splashback radius, manifests itself 
as a sharp density drop in the halo outskirts, at a location that depends upon the mass accretion rate. We 
present calibrations of Rsp and the enclosed mass, Mjp, as a function of mass accretion rate and peak height. 
We find that Rjp is in the range a; 0.8 - IRioom for rapidly accreting halos and is a; 1.5R200m for slowly accreting 
halos. Thus, halos and their environmental effects can extend well beyond the conventionally defined “virial” 
radius. We show that and Rjp evolve relatively strongly compared to other commonly used definitions. In 
particular, Mjp evolves significantly even for the smallest dwarf-sized halos at z = 0. We also contrast Mjp 
with the mass enclosed within four scale radii of the halo density profile, M< 4 r^, which characterizes the inner 
halo. During the early stages of halo assembly, and M< 4 r^ evolve similarly, but in the late stages M< 4 rj 
stops increasing while M^p continues to grow significantly. This illustrates that halos at low z can have “quiet” 
interiors while continuing to accrete mass in their outskirts. We discuss potential observational estimates of the 
splashback radius and show that it may already have been detected in galaxy clusters. 
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1. INTRODUCTION 

In the standard paradigm for structure formation, galaxies 
form through the dissipative condensation of baryons at the 
centers of bound clumps of dark matter, called halos, that 
form within nodes of the cosmic web (Rees & Ostriker 1977; 
White & Rees 1978; Lall & Efstathiou 1980; Blumenthal et al. 
1984). The growth rate of halos controls the rate at which 
baryons are accreted, and complex astrophysical processes 
determine how efficiently these baryons are converted to stars 
(see, e.g., Kauffmann et al. 1993). Therefore, scaling rela¬ 
tions between the luminosity or stellar mass of galaxies and 
the mass of their dark matter halos can provide physical in¬ 
sights into the process of galaxy formation (see, e.g., Kravtsov 
et al. 2004; Tasitsiomi et al. 2004; Vale & Ostriker 2004; Con¬ 
roy et al. 2005; Mandelbaum et al. 2006; Com'oy & Wech- 
sler 2009; More et al. 2009, 2011, 2014; Cacciato et al. 2009, 
2013; Moster et al. 2010; Behroozi et al. 2010, 2013; Zehavi 
et al. 201 1; Leauthaud et al. 2012; Kravtsov 2013; McCracken 
et al. 2014; Kravtsov et al. 2014). Such inference, however, 
requires a basic understanding of the rate of growth of both 
the halos and stellar components of galaxies. 

The halo boundary within which halo mass is measured is 
usually defined as a radius, R/s, of a sphere enclosing a certain 
density contrast A with respect to a chosen reference density 

Pref, 

Mk = ^TrR^Apref . (1) 

Various values of A are commonly used, some motivated 
by the analytic solution for the “virialization overdensity” in 
models of the collapse of a top-hat spherical density perturba¬ 
tion, others motivated by the extent to which observations can 


reliably measure masses (such as Risooc)- Common choices 
for the reference density include the mean matter and critical 
densities of the universe, both of which evolve with redshift. * 

However, the collapse of realistic density peaks in CDM 
models is considerably more complex than is envisioned in 
the spherical top-hat collapse model. Lirst, the peaks are 
not spherical and their collapse rate depends on tidal forces 
which, in turn, depend on the shape of the peak (e.g.. Bond & 
Myers 1996; Dalai et al. 2008a). Second, although the density 
within the peaks is nearly flat at the center (i.e., resembles a 
top-hat profile), it systematically decreases (on average) with 
increasing distance from the peak center (e.g., Bardeen et al. 
1986; Dalai et al. 2010). The collapse of such peaks is thus 
extended in time and is not characterized by a well-defined 
“virialization” epoch. Moreover, the successive collapse of 
density shells results in an extended and smooth density pro¬ 
file. Linally, the peaks have substmcture on smaller scales, 
which results in the collapse and merging of smaller peaks 
during the collapse of a given global peak. Such mergers are 
accompanied by non-linear interactions and the redistribution 
of mass from small to large radii (Valluri et al. 2007). In par¬ 
ticular, in major mergers a significant fraction of the progeni¬ 
tor material ends up at radii beyond the commonly used halo 
boundary (Kazantzidis et al. 2006), meaning that mass is not 
additive in halo mergers when standard mass definitions are 
used. 

A further complication is that Ry, for a given choice of A at 
some epoch zi may not be sufficiently large to enclose all of 
the mass accreted by a halo prior to that epoch. In this case, 

* Although in the case of ACDM model, the evolution of = pdz) satu¬ 
rates at z < 0 as the energy density starts to be dominated by the cosmological 
constant. 
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the change of mass between zi and some later epoch zi < Zi 
will include both the new mass accreted between zi and zi 
and the mass that was accreted at z > zi, but was located 
at r > Thus, in the absence of any actual accretion, 

grows due to the decreasing reference density, but the re¬ 
sulting change in mass cannot be interpreted as physical mass 
growth, change in the halo potential, or as an indication of 
ongoing accretion. 

A number of recent studies (see, e.g., Prada et al. 2006; 
Diemand et al. 2007; Cuesta et al. 2008; Diemer et al. 2013; 
Zemp 2014; Diemer & Kravtsov 2014) have argued that, on 
average, halos of mass < IO'^/z^'Mq accrete little new mass at 
low z at radii < /?a if A > 200 - 300 times the critical density 
of the universe. The increase in their mass Ma is largely due to 
the change of their boundary /?a in response to a decreasing 
reference density, not due to the physical accretion of new 
matter within /?a- In Diemer et al. (2013), we quantified the 
amount of halo mass growth that can be attributed to such 
pseudo-evolution, and found that it contributes significantly to 
the overall mass growth of low-mass halos since z ~ 1. Given 
that standard mass definitions are subject to pseudo-evolution, 
it raises an important question: do galaxy-sized halos accrete 
new matter at late times and if so, at what rate? 

In this paper, we carefully consider particular choices for 
the halo radius and mass definition. We argue that the most 
natural and physical halo boundary can be identified with the 
radius at which newly accreted matter is reaching its first 
orbital apocenter after its initial turnaround. This splash¬ 
back radius corresponds to the outer caustic in the spherical 
models of secondary collapse (Fillmore & Goldreich 1984; 
Bertschinger 1985; Adhikari et al. 2014). The splashback 
radius also physically separates the region where matter is 
infalling for the first time and the region occupied by mat¬ 
ter that has orbited through the central halo region at least 
once. The mass within the splashback radius is thus guaran¬ 
teed to include all of the mass that was accreted by a given 
redshift z. As shown recently (Diemer & Kravtsov 2014, see 
also Adhikari et al. 2014), the halo density profile exhibits a 
sharp steepening of its slope around the splashback radius, 
in correspondence to the density jump expected at this ra¬ 
dius in analytic collapse models (Fillmore & Goldreich 1984; 
Bertschinger 1985; Adhikari et al. 2014). We compare the 
evolution of the splashback radius, R^p, and the correspond¬ 
ing mass, Msp s M{< R^p), with standard mass definitions 
such as Mvir, M 200 C, and M 2 oom- Furthermore, we contrast 
Msp with the mass within a fixed multiple of the halo scale ra¬ 
dius, M< 4 r,. This mass definition characterizes the evolution 
of the inner regions of halos and is manifestly unaffected by 
pseudo-evolution. 

The paper is organized as follows. In Section 2 we discuss 
analytical models of peak collapse and motivate our choice 
of the splashback radius as the physical boundary of a halo. 
In Section 3 we compare the splashback radius and mass to 
several commonly used mass definitions, as well as the mass 
within a fixed multiple of the halo scale radius. In Section 4 
we discuss several implications of the existence of the splash¬ 
back radius, as well as its potential observational signatures 
and possible detections in existing observations. We summa¬ 
rize our results and conclusions in Section 5. In Appendix A 
we present a detailed analysis of various contributions to the 
growth of halo masses. 

Throughout the paper, we denote the mean matter density 
of the universe at the redshift of analysis as and the critical 
density as pc- Mass definitions using a constant overdensity 


A relative to pm or pc are denoted as MAm = M(< R^m), e.g. 
M 2 oom, or Mac = M(< /?ac), e.g. M 2 ooc- Mvir and R„i, denote 
masses and radii defined using the redshift-dependent “virial” 
contrast A(z), computed using the approximation of Bryan & 
Norman (1998). 

2. THE COLLAPSE OL DENSITY PEAKS AND HALO 
MASS 

In this section, we consider the most important processes 
that occur during the collapse of density peaks in hierarchical 
structure formation models, and discuss how these processes 
relate to definitions of the halo boundary and mass. We first 
establish this connection in the context of simplified analytical 
models, and then show that the predictions of these models 
manifest themselves in simulated CDM halos. 

2.1. The Halo Mass and Boundary in the Spherical 
Collapse Model 

In order to elucidate the connection between the collapse of 
a halo and its boundary, we wish to consider a simple, analyt¬ 
ical model. The simplest such model describes the collapse of 
a top-hat perturbation. However, as laid out in the introduc¬ 
tion, peaks in a Gaussian density field have density profiles 
that decrease with radius, significantly different from a con¬ 
stant density spherical top-hat. While a top-hat perturbation 
collapses at a well-defined moment in time, the radial shells 
associated with realistic peaks collapse at different times, re¬ 
sulting in a halo formation history that is extended in time. 
Thus, we expect better guidance from collapse models that 
describe the extended collapse of matter. 

One such model considers the secondary infall of matter 
onto a pre-existing overdensity in an Qn, = 1 universe (Fill¬ 
more & Goldreich 1984; Bertschinger 1985). In this spherical 
collapse model, all mass shells are bound to the pre-existing 
overdensity. Each shell will initially expand with the Hubble 
flow, decelerate, eventually turn around and start contracting. 
At some point, the shell will cross previously collapsed shells 
that are now oscillating in the perturbation potential, thereby 
entering the multi-stream region of the halo. The matter of 
the shell will eventually pass through the pericenter of its or¬ 
bit in the inner region of the perturbation and expand to the 
apocenter of its first orbit. Each successive shell collapses 
onto a deeper potential well than the preceding shell, and thus 
acquires a higher energy and a larger orbit apocenter. 

In this picture, material piles up near the apocenter due to 
its small radial velocity in this region of the orbit, creating a 
density enhancement or caustic which is extremely sharp in 
the case of spherical symmetry (see, for example, the detailed 
discussion in Mohayaee & Shandarin 2006). This caustic oc¬ 
curs at radii r a; 0.1 - 0.4 times the turnaround radius of the 
material at apocenter (see, e.g., Eig. 4 of Vogelsberger et al. 
201 1), depending on the slope of the density profile of the ini¬ 
tial perturbation which determines the mass accretion rate of 
the collapsing halo. Even in the case of ellipsoidal collapse, 
the caustic region is marked by a sharp jump in the density 
profile (Adhikari et al. 2014). 

We note that the outermost caustic corresponds to the apoc¬ 
enter of matter on its first orbit, i.e. the splashback radius 
of the newly accreted matter, which we denote as /?sp- The 
splashback radius cleanly separates the multi-stream region 
of the perturbation at r < R^p from the infall region at r > /?sp, 
where successive shells of matter have not yet crossed. The 
spherical collapse model thus motivates R^p as a natural def¬ 
inition of the halo boundary, and the halo mass as the mass 
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within this radius, Msp s M(< Rsp). By dehnition, the in¬ 
crease in Msp between two epochs Zi and Zf is entirely due to 
the mass shells that have entered the multi-stream region in 
the interval Az - Zi- Zt- Therefore, dM^pjdt is the true halo 
mass growth rate in this model. In contrast, if we had chosen 
a smaller radius < ^sp as the halo boundary, the halo mass 
growth rate would be due to both the accretion of new mat¬ 
ter during Az and matter previously accreted at z > Zi, as we 
discuss in detail in Appendix A. 

2.2. The Mass and Boundary of Realistic CDM Halos 

Although the collapse of realistic CDM halos is consid¬ 
erably more complicated than the collapse of a single peak 
in the secondary infall model, we can still use this model to 
guide our choices of the halo boundary and mass definitions. 
As shown by Diemer & Kravtsov (2014, see also Adhikari 
et al. 2014; Lau et al. 2014), halos that accrete mass at a suf¬ 
ficiently high rate do exhibit a sharp steepening of the den¬ 
sity profile in the outer regions which is due to the caustic 
formed by recently accreted matter. The radius at which the 
profile achieves its steepest slope depends on the halo accre¬ 
tion rate and varies from Rjp ~ O.SRaoom for fast-accreting 
halos to Rsp ~ l -5R200m for slowly accreting halos (Diemer 
& Kravtsov 2014). Adhikari et al. (2014) have conhrmed this 
result and showed that the location of the steepest slope can 
be reproduced using the radius of the outermost caustic in the 
simple model of spherical collapse discussed in Section 2.1. 
In particular, they demonstrated that Rsp depends on both the 
mass growth rate and cosmological parameters such as flni(z). 

Figures 1 and 2 illustrate the correspondence between den¬ 
sity profiles and Rjp using the example of two individual, 
cluster-sized halos with similar masses but very different mass 
accretion rates, representative of the slow and fast accreting 
sub-populations. We operationally dehne the mass accretion 
rate the same way as in Diemer & Kravtsov (2014), 

r= Alog(M™-)/Alog(fl). (2) 

Figure 1 shows the density distribution in a slice through the 
halo center, while Figure 2 shows the spherically averaged 
density profiles and their logarithmic slope. Both hgures con¬ 
trast Rsp (dashed lines) with Raoom (dot-dashed lines) and the 
“virial” radius Ryir (solid lines). The Rsp radii shown in the hg¬ 
ures were predicted using the median relation given by Equa¬ 
tion (5) below and the F of the specihc halos as determined 
from the halo catalogs. Figure 1 clearly demonstrates that the 
density helds exhibit a sharp jump at R ^ Rsp, and that this 
radius occurs at a smaller multiple of R 200 m for the faster ac¬ 
creting halo. We conhrm these impressions by considering 
the spherically averaged density prohles of the same halos in 
Figure 2 which highlights how steep the density profile can 
get around Rsp (a logarithmic slope of -7). 

The correlation between Rsp/Rzoom and the accretion rate 
mimics the correlation between the ratio of the last caustic 
and turnaround radii and the slope of the initial perturbation 
profile in the spherical collapse model (see Vogelsberger et al. 
2011, and Section 2.1). By further analogy with the spherical 
collapse model, the splashback radius of CDM halos should 
include approximately all of the mass ever accreted by a halo. 
Thus, changes in Msp should always correspond to the current 
accretion of new mass, implying that Msp is largely unaffected 
by pseudo-evolution. In Sections 3.1 and 4.6 we discuss how 
Rsp can be measured in cosmological simulations and obser¬ 
vations, and how Msp and its evolution relate to conventional 
spherical overdensity masses. 


3. HALO RADIUS AND MASS EVOLUTION 

In the previous section we argued that the most natural and 
physical definition of the mass associated with a density peak 
is the mass enclosed within the radius of the outermost caus¬ 
tic, Rsp, which we now consider in detail. Eollowing the dis¬ 
cussion of Rsp and Msp and their calibrations using cosmolog¬ 
ical simulations, we also consider a halo boundary defined as 
a constant multiple of the scale radius, R = 4rs, and the cor¬ 
responding mass, M< 4 r^, which characterizes the mass in the 
inner regions of halos. We then compare the evolution of halo 
mass and radius in these definitions to the commonly used 
“virial” mass definitions. 


3.1. The Splashback Mass, M^p 

As shown by Diemer & Kravtsov (2014), Rsp can be mea¬ 
sured from the halo density profile as the radius where the 
density profile steepens sharply beyond what is expected from 
the NEW and Einasto predictions. We use the same simu¬ 
lation suite and fit the median density profiles of halo sam¬ 
ples with a range of different masses, accretion rates, and red- 
shifts using the fitting function in Equation (4) of Diemer & 
Kravtsov (2014), which we reproduce here for completeness. 


Pif} ~/trans PEinasto + Pouter 


PEina.sto ~ Ps CXp 


ftians ~ 


1 + 1 - 


-r//3 


Pouter “ Prr 


5R 


200m 


+ 1 


(3) 


Eor consistency, we fix some of the parameters in the fitting 
function in all fits (regardless of whether the samples were 
selected by mass or accretion rate), namely /? = 6, y = 4. 
Similar to Diemer & Kravtsov (2014), we also fix a according 
to the relation with v as calibrated by Gao et al. (2008), 

Q'(v) = 0.155 + 0.0095(4) 


The other parameters, namely rs, rt, he, and Se, are determined 
from a least-squares fit. 

Diemer & Kravtsov (2014) showed that, at z = 0, the 
turnover radius, rt, at which the density profiles steepens, de¬ 
pends on the mass accretion rate F. In Figure 3, we extend this 
analysis to higher redshifts and use Rsp, dehned as the radius 
where median profile of halos reaches the steepest slope, in¬ 
stead of rt. The redshift intervals over which F are measured 
are the same as the redshifts listed in Figure 3, i.e. for z = 0, 
F is measured between z = 0.5 and z = 0, for z = 0.5 between 
z = 1 and z = 0.5, and for z = 4 between z = 6 and z = 4. 
The choice of the redshift intervals defining F is somewhat 
arbitrary, but corresponds reasonably closely to the expected 
crossing time through the full extent of the halo, 2R (for ex¬ 
ample, at z = 0, Az corresponds to about 5 Gyr whereas the 
crossing time is about 4 Gyr). 

We bin halos both by v and by F, and only use halo samples 
with V > 1 (M > 3 X 10'2/!-'Mo atz = 0, M > lO" h-^M^ at 
z = 1) for this analysis, as the density jump associated with the 
splashback radius is difficult to measure robustly from spher¬ 
ically averaged prohles in halos with low v and low F. This 
issue is apparent in Figure 10 in Diemer & Kravtsov 2014: for 
prohles with low v and low F, the 2-halo term begins to dom¬ 
inate at radii smaller than Rsp, thus concealing the steepening 
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Figure 1. Projected density in a slice of thickness 0.15/?200ra through the center of two halos with low (left, F = 0.8) and high (right, F = 2.7) mass accretion 
rates. The halos have similai* masses, Mvir = 1-1 x 10^"^ and 1.8 x 10^"^ h~^M q at z = 0. The white lines show 7?vir (solid), i? 200 m (dot-dashed), f?sp (dashed) 
and i?infaii (dotted; see §3.1 for a detailed description of these radii). /?sp and f?infaii were calculated using the calibrations presented in Section 3.1 rather than 
the density profiles of the individual halos shown. Halos with a low mass accretion rate exhibit a caustic at a radius significantly larger than /? 200 m> whereas 
fast-accreting halos have i?sp ~ ^ 200 m (at z = 0). The visualizations were created using the algorithm of Kaehler et al. (2012). 



T /^200m 


/ -R200ni 


Figure 2. Spherically averaged density profiles (top panels) and their logarithmic slope (bottom panels) of the two halos shown in Figure 1. The slopes were 
computed using a profile smoothed with the fourth-order Savitzky & Golay (1964) filter over the 15 nearest bins. The steepening around f?sp is very pronounced 
in both profiles, but the profile of the faster accreting halo reaches a steeper slope and at a smaller radius. The vertical lines in the bottom panels mark the same 
radii shown in Figure 1 using the same line types, i.e. ^sp. and /?infaii (defined as the radius where the mean radial velocity profile of Vr reaches minimum) 
from left to right. For the slower accreting halo (left), the estimate of Equation 5 slightly underestimates the true i?sp- This disagreement is not surprising since 
the /?sp of individual halos are expected to scatter around the median relation. 
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Figure 3. Splashback radius, 7?sp (top panel), and the mass within this radius, 
Msp (bottom panel), as a function of the halo mass accretion rate, F. Darker 
points correspond to halo samples with higher peak height, v. The halos 
were binned in v-bins of width 0.5, starting at v = 1. Samples with v < 1 
were omitted as R^p is hard to measure for their profiles (see the discussion in 
Section 3.1). The figure demonstrates that/?sp/^ 200 m and Msp/M 200 m depend 
on r and z, but do not show a strong dependence on v at fixed accretion rate 
and redshift. For halos in our ACDM cosmology these dependencies can 
be approximated by Equations (5) and (6), shown with a solid line for each 
redshift. 

in the density profile. This does not mean that low-mass halos 
do not exhibit a steepening in their density profile; however, 
they are strongly influenced by their environment, making it 
difficult to discern the location of /?sp- 
Figure 3 demonstrates that, at fixed F, Rsp/^ 200 m does not 
depend on v, but does depend on z. In particular, Adhikari 
et al. (2014) showed that the overdensity associated with 
the splashback radius depends on where Qm(2) = 

Pm{z)lpmtiz)- Thus, we parameterize the dependence of 
f^sp/f^ 200 m on F and z with the fitting function 

= 0.54 [1 -H 0.53Q„(z)] (l -H 1.36e-r/^°4) , (5) 

^200m ^ 



0 1 2 3 4 5 

^^200m 

Figure 4. Splashback radius, i?sp (top panel), and the corresponding enclosed 
mass, Msp (bottom panel), relative to ^ 200 m and M 200 m as a function of peak 
height, V 200 m = IO'i^ioom)ID{z). The dotted vertical lines indicate v = 
1, the peak height below which the relation was not directly calibrated (see 
Figure 3). The dependence of ^sp/^ 200 m aiises because halos with higher v 
have, on average, a higher mass accretion rate, and thus a smaller ^sp/^ 200 m 
(Figure 3). The dashed lines show the relations specified in Equations (7) 
and (8). The figure shows that for the rare, massive halos that accrete mass 
at a fast rate the splashback radius is close to /? 200 m, while for low-v halos it 
extends to a considerably larger radius of R^p ~ 1.5i?200m- 

shown with solid lines in Figure 3. Given this function, we 
could now compute the median Msp from Rjp by assuming a 
particular form of the density profile. However, we get a more 
accurate fit by directly calibrating the median ratio of Msp and 
M 2 oom using the simulated density profiles, 

-2^ = 0.59 [1 -H 0.35Q.^(z)] (l -H 0.92e-r/4.54\ 

M200m ' ' 

These formulae were calibrated using a flat ACDM cosmol¬ 
ogy where Qm = 1 - Qa = 0.27, h - 100 /Hq - 0.7, 
CTg = 0.82, tis = 0.95. These calibrations can be used to com¬ 
pute the dependence of the overdensity. A, of halos on F and 
Qtn(z). Since we do not find a large dependence of Rsp/f^ 200 m 
on V, we will extrapolate the relations calibrated above even 
for V < 1 halos in the subsequent sections while discussing 
our results. We also note that Adhikari et al. (2014) have pre¬ 
sented a calibration for A as a function of the instantaneous 
mass accretion rate, s and Dni(z). Once the differences be¬ 
tween s and F are accounted for, our calibration is largely 
consistent with theirs. 

In observations, however, the accretion rate or the exact 
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density profile of a halo are not readily available. Thus, we 
also quantify the dependence of R^p and Msp on the conven¬ 
tionally defined, observable Ma, or rather peak height, v s 
6cla-{M/^)ID{z), in Figure 4.^ This dependence arises because 
halos of higher peak height exhibit, on average, higher accre¬ 
tion rates (see e.g.. Figure 8 in Diemer & Kravtsov 2014). 
In order to translate Equations (5) and (6) into functions of v 
rather than F, we use the model of Zhao et al. (2009) to cal¬ 
culate halo mass growth histories. For each redshift along an 
accretion history we compute the accretion rate F across the 
same redshift intervals as in Diemer & Kravtsov (2014) and 
calculate R^p and M^p using Equations (5) and (6). Eigure 4 
shows the results as a function of peak height. The relations 
are more or less independent of redshift, and well fitted by the 
approximations 


R.p 

= 0.81 

(l -b 0.97e"'’''2 '^‘*) 

(7) 

RlOOm 

-^sp 

= 0.82 1 

(l H-0.63e^’'/^-^2) , 

(8) 

47200m 


shown with dashed lines in Eigure 4. Thus, in the con¬ 
cordance cosmological model, the dependence of f?sp/^ 200 m^ 
4^sp/4^200m or the overdensity of halos on Qm)^) is approxi¬ 
mately cancelled by the dependence of F(v) on redshift based 
on the median mass accretion histories of halos. Eor rare, 
massive halos that accrete mass at a fast rate (on average), the 
splashback radius is close to /? 200 iin while for low-v halos it 
extends to a considerably larger radius of R^p ~ 1.5/?200m- 

The R^p and Msp calibrations presented in Equations (7) and 
(8) were obtained using the median profiles of halos of a given 
F or V. In observations, however, stacking would result in an 
average of the density profile, not the median R^p and Mjp. 
We have checked that the v-R^p and v-Msp relations obtained 
from averaged profiles are almost identical to Equations (7) 
and (8). 

Given that real CDM halos are not spherical, the density 
jump associated with the splashback radius occurs at different 
radii in different directions from the halo center. The corre¬ 
sponding feature in the spherically averaged density profile is 
thus not a sharp jump but rather a steepening of the profile that 
spans a range of radii. The finite radial extent of the steepen¬ 
ing creates a certain ambiguity in the choice of the splashback 
radius definition. Eor instance, an alternative definition of the 
splashback radius could be the radius where the average ra¬ 
dial velocity in a shell is most negative, /^infaii. This radius is 
more likely to include most of the accreted mass, although the 
majority of the mass between R^p and /?infaii is infalling for the 
first time. We find that /?infaii ~ 1 -T/^sp and Minfaii ~ 1.2Msp 
at all redshifts and halo masses. Thus, the fitting formulae in 
Equations (7) and (8) can easily be modified to return f?infaii 
and Minfaii. 

Eigure 1 above shows both R^p and /?infaii ~ 1 •4/?sp (dotted 
line). In contrast to R^p, /^infaii clearly extends into the fila¬ 
mentary regions and is not associated with the collapsed halo 

^ The peak height is defined using the M 2 oom mass, i.e. v = V 2 oom = 
Sc /(T(M2oom)/D(z), where Sc is the critical threshold for collapse, cr^ is the 
variance of initial density fluctuations when smoothed with a top-hat filter 
with a size corresponding to the Lagrangian radius of mass M 2 oom, and D(z) is 
the growth factor. For the color scale of Figure 3, we use Vyir for compatibility 
with Diemer & Kravtsov (2014). Note, however, that the difference between 
Vvir and V 200 m is < 5% for all masses and redshifts. See Diemer & Kravtsov 
(2014) for the exact definition of cr{M). 


matter. Eigure 2 confirms this impression, as /^infaii does not 
correspond to a particular feature in the density profiles. Thus, 
the splashback radius definition based on the steepest density 
profile slope is preferable, and will be used for the remainder 
of this paper. 

3.2. The Inner Mass, 

The results presented in Section 3.1 demonstrate that the 
splashback radius of halos is quite large. Eor some purposes, 
it may be instructive to consider the mass evolution in the 
inner regions of halos. In principle, one could character¬ 
ize the inner regions simply by using Ma with a high value 
of A. However, given that any spherical overdensity radius 
Rh < ^sp is subject to pseudo-evolution, we would prefer to 
define the inner mass using a radius that is not tied to any 
cosmological reference density. 

Zemp (2014) discuss several alternative halo mass and ra¬ 
dius definitions that are not subject to pseudo-evolution. In 
particular, it was argued that the mass within a radius that en¬ 
closes a fixed physical density independent of redshift, such 
as 200prn(2 = 0), could be used as an alternative measure 
of mass. Such a definition, however, has a number of draw¬ 
backs. Eirst, if the threshold density is chosen to be too low, 
the corresponding radius will be much larger than the virial- 
ized region of halos at high z. If the density is chosen to be 
too high, the enclosed mass will correspond only to the inner 
region of the halo. Most importantly, the mass and radius de¬ 
fined in this way do not track the physical growth of a halo in 
its fast accretion regime, where the halo profile does exhibit 
a well-defined characteristic scale close to /? 200 m( 2 ) (Diemer 
& Kravtsov 2014, see also Eigure 5). Zemp (2014) has also 
discussed the possibility of using the scale radius (the radius 
where the density profile has a logarithmic slope of -2) as the 
halo boundary, and define the mass as M(< rj). 

Here we use a halo radius and mass definition where R - 
4rs and M< 4 r^ = M(< 4rs). The multiple of 4 in this def¬ 
inition is motivated by the fact that the concentration, ca = 
T^a/g, is approximately equal to four as long as the halo is 
in the fast accretion regime (Zhao et al. 2003, 2009).^ In 
this regime, where the contribution of pseudo-evolution to the 
mass growth is relatively small, M< 4 r^ approximately tracks 
Ma. Subsequently, as the mass growth and physical evolution 
of the inner region of the halo profile slow down, the scale 
radius approaches a constant (Bullock et al. 2001) and M< 4 r^ 
tracks the actual evolution of the inner halo mass due to real 
profile changes, unaffected by pseudo-evolution. Assuming 
an NFW density profile, M< 4 r^ is given by 

y«(4) 

M<4a=Ma-^. (9) 

R{ca) 

where jj(x) - ln(l + x) - x/{l + x) and the mass and concen¬ 
tration could correspond to any of the commonly used density 
contrast choices. 

Eor example, for a Milky-Way sized halo of mass Myir = 
10'^ /z^'Mq (f?vir = 207 /j^'kpc) and a typical concentration at 
that mass, Cyir ~ 9 (e.g., Zhao et al. 2009; Diemer & Kravtsov 
2015), 4rs 92 /i^'kpc, and M< 4 r^ 5.8 x lO" /z^'Mq. Eor 

comparison, the median R^^p for such a halo is 358 h^^kpc « 
511 kpc and the median Msp is 1.4 x 10^^ about 2.4 

^ The concentration in the fast accretion regime has a residual dependence 
on halo mass and redshift and thus varies between Cmin » 3 and Cmin » 4 
(Diemer & Kravtsov 2015). However, for the sake of simplicity, we choose a 
fixed value of Cn,i„ = 4. 
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Figure 5. Median mass growth histories (top) and halo boundary growth rates (bottom) of halos of different masses Ma at z = 0 (solid lines), as well as the 
mass within the splashback radius (dashed lines), calculated as described in the text. In the left column we use Ma = M 200 C. in the middle column Ma = Myir, 
while in the right column Ma = M 200 m- In the fast mass growth regime, Msp « Myii-, while in the slow mass growth regime Msp grows faster than Myir- The 
transition between these two regimes happens earlier for smaller halos. The differences between M^p and other mass definitions are small at high redshift, but 
increase with decreasing redshift. For low mass halos, the high contrast mass definitions such as M 200 C can show differences of up to 50% at z = 0. For cluster 
mass halos Msp % M 200 m- The differences in the halo boundaries are much larger than those in the halo masses. 


times larger than M< 4 r^. We caution that these values are me¬ 
dians, and there is large scatter both in the concentration-mass 
relation and in R^p at fixed mass. 

3.3. The Evolution of Halo Radii and Masses in Different 
Definitions 

We now contrast the redshift evolution of R^p and M^p with 
the evolution of the boundary and mass for the definitions 
Ma, using some common choices of A. Once again, 
we use the model of Zhao et al. (2009) to calculate concentra¬ 
tions and halo mass growth histories. The concentrations are 
used to convert between the different choices of A, while the 
halo mass growth histories are also used to derive T (and thus 
Rsp and Msp using the fitting formulae presented in Section 
3.1). 

The upper panels of Figure 5 show the evolution of Msp and 
the traditional definitions M 200 C. A^vir and M 2 oom- At high z 
(i.e., in the fast mass growth regime), Msp ^ M 200 C ~ ^vir ~ 
A^ 200 m- At low z, Msp evolves faster than Myir for halos of 
all masses. In particular, the figure demonstrates that even 
dwarf and Milky Way-sized halos do accrete new mass at low 
redshifts. In the lower panels, we compare the evolution of 
the corresponding radii. While Rsp can be significantly larger 
than R 200 m (see also Figures 3 and 4), only a relatively small 
fraction of the total halo mass resides at those radii, reducing 


Msp/M 2 oom compared to Rsp/f^ 200 m- 

The differences in the growth rates of Msp, M 200 C. Mvk, and 
Mioom arise because the mass accreted within Rsp at late times 
is distributed with an approximately isothermal p oc pro¬ 
file and contributes significantly only at large radii, r > Ryi^ 
(see Section 4.1). This highlights an important point; at low 
redshift, the halo mass distribution in the inner regions may 
be relatively stable and evolve slowly, while the outer regions 
may evolve fast. 

It is therefore also interesting to contrast the evolution of 
Msp with M< 4 r^ which characterizes the mass distribution in 
the inner regions of halos at low z (see Section 3.2). This 
comparison is shown in Figure 6, both for the masses (top 
panel) and radii (bottom panel). The most massive halos are 
largely still in the fast accretion regime today, and the evolu¬ 
tion of Msp and M< 4 r, are quite similar. Low-mass halos, on 
the other hand, are in the slow mass growth regime and their 
Msp and M< 4 r^ evolve quite differently. For example, for ha¬ 
los of Mvir(z = 0) = 10* /i^'Mq, M< 4 r^ evolves significantly 
slower than Msp(z) at z < 3.5. At z < 1, M< 4 r^ is approxi¬ 
mately constant, while Msp for these halos changes by « 30% 
between z = 1 and z = 0. Finally, we compare the halo mass 
growth rates for the Msp, M 200 C and M< 4 r^ definitions in Fig¬ 
ure 7. The growth rate of M< 4 r, is approximately an order of 
magnitude lower than that of Msp at z < 1 for low-mass halos. 













8 


More, Diemer & Kravtsov 



2 


Figure 6. Median mass growth (top) and halo boundary growth histories 
(bottom) of halos of different virial masses at z = 0 (solid lines) for the splash¬ 
back definition, as well as for the mass within 4 scale radii, (dashed 

lines), calculated using the model of Zhao et al. (2009). In the fast growth 
regime, « Mvir ~ Msp, while in the slow growth regime M< 4 rj. ap¬ 

proaches a constant value. The transition between these two regimes happens 
earlier for smaller halos. The slight decrease in at low z appears to be 

a small artefact of the Zhao et al. (2009) model and is not present when we 
plot a similar evolution using halos from cosmological simulations. 

In summary, Figures 5-7 show that during the fast accre¬ 
tion regime Myir, M 2 oom. and M< 4 r^ are approximately 
equivalent. In the slow mass growth regime, however, the 
mass within the inner radii (r < r^r^) for halos with < 
10'^ hr^MQ stops growing, while Mjp keeps growing even 
faster than Myir. 

4. DISCUSSION 

In the previous sections, we have argued that a natural def¬ 
inition of the halo boundary is the splashback radius, which 
corresponds to the apocenter of matter on its first orbit after 
accretion. Although the distribution of the apocenters of re¬ 
cently accreted material is not spherical as in the idealized 
models of secondary accretion that motivate the concept, real 
CDM halos do exhibit a steepening of their radial density 
profiles over a relatively natTow range of radii. We dehne 
the splashback radius as the radius of the steepest slope of 
the density profile and calibrated this radius and the mass en¬ 
closed within it for the concordance cosmological model. Us¬ 
ing the splashback radius as a halo boundary has a number of 
implications. 



z 

Figure 7. Comparison of the growth rates of the splashback mass, Mgp, 
M< 4 ^^, and M 200 c- At low redshifts the difference between the rates based 
on M 200 C (^<4rs) ^sp can be as large as a factor of 2 (10) for low-mass 
halos. 

For instance, the fact that f?sp/^ 200 m Msp/M 2 oom are in¬ 
dependent of redshift at a fixed v (Figure 4) may have impli¬ 
cations for the universality of the halo mass function. The 
mass function is found to be most universal when masses are 
defined with respect to the mean density of the universe, such 
as M 2 oom. and less universal when masses are dehned with 
respect to the critical density. From general considerations, 
one may expect peak collapse process and the associated halo 
mass function for Mjp to be approximately universal. The uni¬ 
versality of Msp/M 2 oom as a function of v would then explain 
why the halo mass function for the M 2 oom definition is ap¬ 
proximately universal, or at least why it is considerably more 
universal than for the M 200 C definition. 

We have shown that, at late times, the evolution of Rsp and 
the enclosed mass, Mjp, is considerably faster than that of 
^ 200 c or /?yir and their corresponding masses. However, 
which characterizes the mass distribution in the inner regions 
of slow accreting halos, evolves very little at z < 1. This dif¬ 
ference implies that a quiescent evolution of the inner regions 
of halos can co-exist with active growth in the outer regions. 
Less massive halos at low redshifts, in particular, do accrete 
mass, even though their inner regions bear scant evidence for 
such growth. We have also shown that for slow accreting ha¬ 
los at low redshifts the splashback radius can be up to a factor 
of two larger than /? 200 c. which is often used to define the halo 
boundary in galaxy formation studies. This raises the possi¬ 
bility that the “zone of influence” of individual halos (where 
satellite-specific environmental effects can be seen) may be 
much larger than is usually assumed. In the remainder of this 
section, we discuss these issues, as well as operational defi¬ 
nitions of /?sp and M< 4 r, in simulations and observations and 
possible observational detections of the splashback radius. 

4.1. Quiet Interiors with Active Outskirts 

Figures 5 and 6 show that M^p grows signihcantly all the 
way to z = 0 for halos of all masses, while the mass in the 
inner regions (characterized by M< 4 ,.J evolves considerably 
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Figure 8. Median mass growth (top) and of cluster-sized halos for mass 
within two NFW scale radius, M 2 t^ (dashed lines) and mass within the splash¬ 
back radius (solid lines) calculated using the model of for mass accretion his¬ 
tory and concentrations Zhao et al. (2009). The figure shows that for 10*'* Mq 
halos M 2 rs grows from z = ltoz = 0bya factor of two, while Msp grows 
by a factor of four. For 10*^ Mq halos which grow at a much faster rate, the 
growth of the two masses over the same redshift interval is comparable. 

slower at z < 1-2, and not at all for galaxy-sized halos. 
The reason for the slow evolution of the inner regions is not 
that the overall accreted mass is small (galaxy-sized halos 
approximately double their Msp between z - I and z = 0; 
see Figure 5), but that the newly accreted matter has a rela¬ 
tively shallow density profile (see Section 4.2 in Lithwick & 
Dalai 201 1). The shallow profile arises because the radial pro¬ 
file of the potential over most of the halo volume is shallow, 
meaning that the velocity of the accreted matter does not vary 
strongly with radius. Hence, the time spent at each radius r 
is 6t(r) ~ rjv oc r. The time averaged mass profile of newly 
accreted mass is thus Macc(< r) ^ 6t{r) oc r/v oc with 
a ~ 1, leading to a density profile which is close to isother¬ 
mal, p oc This profile is shallower than the overall NFW- 
like profile of CDM halos, which implies that the previously 
accreted matter dominates in the inner regions of halos while 
newly accreted mass contributes significantly in the outer re¬ 
gions near the splashback radius. 

This highlights the possibility that systems that appear qui¬ 
escent in their interior regions may actually still actively grow, 
particularly when much of the mass growth is due to the accre¬ 
tion of diffuse mass and small halos rather than due to major 
mergers. The orbit of this material may take it to the interior 
regions, but it is predominantly deposited in the outskirts. Let 
us consider some implications for both massive, cluster-sized 
halos as well as galaxy-sized halos. 

For cluster-sized halos, relaxed “cool-core” inner regions 
do not preclude the possibility of active accretion in the outer 
regions. Abell 133 (further discussed in Section 4.5) may be 
an example of such a system, as it combines the characteris¬ 
tics of a cool core system (albeit a disturbed one) and a sys¬ 
tem undergoing rapid mass accretion. We indicate this pos¬ 
sibility in Figure 8 which contrasts the average evolution of 
mass (calculated using the Zhao et al. 2009 model) within 



Figure 9. The growth in the stellar mass of central galaxies along the progen¬ 
itor history of halos in dilferent bins of present-day halo mass. The shaded 
regions show constraints inferred by Behroozi et al. (2013) using detailed 
models which follow galaxy populations in dark matter halo merger trees 
from numerical simulations. We convert the - Mvii- relations they infer at 
different redshifts to M* - M<4rs ( 2 ), and show the redshift evolution of M* 
at fixed M< 4 rs using solid lines. For clarity, we shift both the solid line and 
the shaded region for the mass bins logM^jJ^ = [11.0,11.2] and [12.0,12.2] 
vertically by 0.1 and 0.2, respectively. 

two NFW scale radii, M 2 r,, and mass within the splashback 
radius, M^p, for cluster-sized halos. We show the mass within 
2r, 2R200c/c200c ~ 0.5R200c ~ ^sooc because this radius is 

often used to measure halo masses in cluster studies. The fig¬ 
ure shows that for 10^"^ /i^'Mq halos M 2 r, grows from z = 1 
to z = 0 by a factor of two, while Msp grows by a factor of 
four. For 10'^ /i^'Mq halos which grow at a much faster rate, 
the growth of the two masses over the same redshift interval is 
comparable, which indicates a significant growth of both the 
inner and outer regions of the halo. We conclude that individ¬ 
ual cluster-sized halos can exhibit very different types of mass 
evolutions: for some halos the difference between the growth 
of M 2 r, and Mjp is more than a factor of two, while for other, 
fast growing halos the difference may be small. 

For galaxy-sized halos, the slow evolution of the inner re¬ 
gions can be exploited to identify progenitor halos across red- 
shifts at z < 1. In particular. Figure 6 shows that M< 4 r^ for 
halos of Mvir < lO'^ Mq evolves very little at z < 1. This 
suggests that changes in galaxy properties at fixed M< 4 r^ can 
be used to deduce the evolution of these properties along their 
progenitor histories. 

As an example, we estimate the stellar mass growth in low- 
mass halos at z < 1. For this purpose, we assume the stel¬ 
lar mass-halo mass relation (SHMR) and its evolution as pa¬ 
rameterized by Behroozi et al. (2013), and convert Mvir(z) to 
M< 4 r,( 2 ) to obtain the M, - relation as a function of red¬ 
shift. We can now read off the star formation history at fixed 
M< 4 r,, shown as solid lines in Figure 9. The shaded areas 
show the observed star formation histories as derived from the 
sophisticated modeling of Behroozi et al. (2013). They con¬ 
vert observational data on the redshift dependence of the stel¬ 
lar mass function and of the specific star formation rate stel¬ 
lar mass relation using halo merger trees, to infer the growth 
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Figure 10. SHMRs for central (solid blue) and satellite (solid red) galaxies, 
as obtained by Watson & Conroy (2013). The dashed red line corresponds 
to their model which posits that the stellar mass in satellites grows by 10% 
compared to centrals after they fall into the host halo. The dashed green cuiwe 
shows the SHMR of centrals, but with halo masses measured at the average 
infall redshift, z = 0.5. Due to pseudo-evolution, the measured halo mass 
decreases at higher z, explaining most of the difference between the central 
and satellite SHMRs. 

histories of stellar mass in halos. On the other hand, our sim¬ 
ple inference at fixed M< 4 r^, only utilizes constraints on the 
stellar mass—halo mass relation at different redshifts, such 
as the ones that can be obtained using abundance matching. 
The stellar mass growth histories agree well with the detailed 
modelling results at low halo masses, and starts to diverge 
once the evolution of M< 4 r, becomes significant. Thus infer¬ 
ences related to the stellar mass growth histories in studies of 
the evolution of the galaxy-halo connection (e.g., Leauthaud 
et al. 2012) at low halo masses can be simplified by casting 
those results as a function of . However, as we will dis¬ 
cuss later, tying the growth histories of galaxies to the baryon 
accretion rate in the central regions is considerably more com¬ 
plicated. 

4.2. The Pseudo-evolution of Halo Mass in the Standard 
Definitions 

As we saw in the previous section, most of matter accreted 
by halos is deposited at the outskirts, while their inner regions 
can be relatively quiet. This finding relates directly to the 
mass accretion rate in different mass definitions, and to the 
question whether this accretion is physical or due to pseudo¬ 
evolution. 

First, we note that the halo masses in the standard defini¬ 
tions are within a factor of ~ 1.5 of Mjp at all z (Figure 5), 
though higher density contrasts (such as Msooc) would lead 
to smaller radii and thus larger differences with Mjp. The 
mass accretion rate of Mjp, on the other hand, differs signifi¬ 
cantly from the spherical overdensity definitions, and can be 
up to three times larger than the change in M 200 C (Figure 7). 
This may seem surprising, since the mass growth in conven¬ 
tional definitions already suffers from pseudo-evolution, i.e. 
a growth in addition to the physical growth within the corre¬ 
sponding halo boundary. 

However, there are two competing differences between 


dMspIdt and dMxooddt. While dM 2 oocldt is larger than the 
physical change of the mass within R 200 C due to pseudo¬ 
evolution, it also misses mass that is added outside R 200 C be¬ 
tween two epochs, i.e. the physical growth in the outskirts 
discussed in Section 4.1. While it is not a priori clear which 
of the two effects should dominate. Figure 7 demonstrates 
that dMspIdt is larger than dM 2 oocldt at all masses and red- 
shifts shown, meaning that physical evolution in the outskirts 
dominates over the pseudo-evolution of M 2 ooc- However, for 
halos which are tmly non-accreting (dMspIdt 0), pseudo¬ 
evolution would still lead to a non-zero dM 2 oocldt. 

Whether the evolution in under- or overestimates the 
true accretion rate, it does not always reflect the physical ad¬ 
dition of matter. The relatively small difference between Ma 
and Msp does not reflect the physical correctness of Ma, but 
rather the shallow mass profile in the outer regions of ha¬ 
los which means that large changes in radius lead to small 
changes in enclosed mass. 

The importance of taking pseudo-evolution into account 
when interpreting changes in halo mass can be illustrated by 
the results of Watson & Conroy (2013) who tested the com¬ 
mon assumption that the same stellar mass-halo mass rela¬ 
tion (SHMR) is valid for both host and satellite halos. They 
showed that the relations are somewhat different and inter¬ 
preted the difference as a 10% growth in the stellar mass of 
satellite galaxies since the time of infall, possibly due to resid¬ 
ual star formation before final quenching (Figure 10). They 
estimated the average time of infall to be around 4 Gyr ago, or 
z ~ 0.5. However, in their analysis (similar to all abundance 
matching studies) the halo masses of centrals are defined at 
z = 0 while those of satellites are defined at their infall epoch 
(z — 0.5). Therefore, some of the difference in the host and 
satellite SHMRs will arise from the pseudo-evolution of mass 
between these redshifts. We estimated this effect and show it 
with the green dot-dashed line in Figure 10. Simply account¬ 
ing for pseudo-evolution explains most of the difference in the 
SHMRs found by Watson & Conroy (2013), and does so bet¬ 
ter than the suggested blanket 10% increase in the stellar mass 
of satellites across all halo masses. This agreement suggests 
that satellite galaxies maintain levels of star formation similar 
to their central counterparts for a substantial period after in¬ 
fall. This increase in stellar mass would depend on a satellite’s 
halo mass, a qualitatively different view from that presented 
by Watson & Com'oy (2013). 

4.3. The Accretion of Mass onto Low-mass Halos 

A number of recent studies (Prada et al. 2006; Diemand 
et al. 2007; Cuesta et al. 2008; Diemer et al. 2013; Zemp 
2014; Diemer & Kravtsov 2014) have pointed out that the in¬ 
ner density profiles of low-mass, galaxy-sized halos evolve 
very little at z < 1. The near-constant density profiles were 
sometimes interpreted as a lack of new mass accretion in these 
halos. As we showed in the previous section, Msp does, on 
average, increase rapidly for halos of all masses, including 
dwarf-sized halos, implying the continued accretion of mass 
in all halos to low z. This newly accreted mass may briefly 
pass through the inner regions of halos, but is primarily de¬ 
posited in the halo outskirts, with no substantial change in the 
inner regions. 

In Section 4.1, we explained why the newly accreted mass 
contributes little to the inner density profile, but this argu¬ 
ment applies only to collisionless dark matter. Baryons, on 
the other hand, experience additional pressure forces and pos¬ 
sibly dissipation and interaction with feedback-driven winds 
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Figure 11. Comparison of the satellite fraction as a function of the peak cir¬ 
cular velocity, Fpeak^ for subhalos identified as density peaks within the virial 
radius (solid line) to those identified within the splashback radius (dashed 
line) of more massive halos. The satellite fraction significantly increases for 
low values of Vpeak- 

(Nelson et al. 2015), so that their radial distribution and fate 
can be quite different (Faucher-Giguere et al. 2011; van de 
Voort et al. 2011; van de Voort & Schaye 2012; Wetzel & Na- 
gai 2014; Feldmann & Mayer 2015). We note that the baryon 
accretion rate in the inner regions of halos could track the 
growth rate of halos (e.g., determined by Mjp) if the trajectory 
of recently accreted material brings it closer to the central den¬ 
sity peak. The accreted baryons could interact with intra-halo 
gas and feedback-driven winds and stay in the inner regions 
of halos. Such effects could conceivably explain the large 
baryon accretion rate and its consistency with f\,dM 2 oomldt, 
as found by Wetzel & Nagai (2014). However, Woods et al. 
(2014) tracked particles in their hydrodynamical simulations 
and showed that the late time star formation in their simu¬ 
lated galaxies is a result of gas cooling from a reservoir which 
was accreted 2-8 Gyr earlier. Nelson et al. (2015), on the 
other hand, find that feedback-driven winds increase the time 
it takes the gas to reach the disk after it crosses Rvir by a factor 
of several, but does not introduce any significant mass depen¬ 
dence in the accretion rate. 

Although the fate of baryons is not completely clear, our re¬ 
sults do show that galaxy-sized systems do have a continuing 
supply of fresh matter to z = 0. 

4.4. A Halo’s “Zone of Influence” and the Extent of the 
One-halo Term 

We have shown that the splashback radius for galaxy-sized 
halos at low z can be more than a factor of two larger than 
Riqoc- For example, for Milky Way-sized halos, the splash¬ 
back radius can be as large as 600 kpc, or two thirds of the 
distance to M31 (see, e.g.. Figure 2 of Diemand & Kuhlen 
2008 and Figure 5 above). The large Rjp may explain the 
presence of dwarf spheroidal galaxies, such as Cetus (e.g., 
Fraternali et al. 2009), at large distances, and highlights that 
the “zone of influence” of halos (i.e., their associated environ¬ 
mental effects) can extend to and beyond 2/?200c- 

This effect was discussed previously in the context of 
“backsplash satellites” which were found in substantial num¬ 
bers beyond the halo radii in the standard definitions (Balogh 
et al. 2000; Mamon et al. 2004; Gill et al. 2005, see also Wet¬ 


zel et al. 2014 for in depth analysis of this issue and more 
recent references). Their occurrence at large radii is not sur¬ 
prising because satellite subhalos effectively trace the distri¬ 
bution of accreted matter. This also means that the radial num¬ 
ber density profiles of subhalos in simulations and galaxies in 
observed clusters may exhibit the steepening at the splash¬ 
back radius (see further discussion in Section 4.5). Consid¬ 
ering the extent of the splashback radius, it is unnecessary to 
invoke special ejection processes to explain the presence of 
subhalos that have orbited their host at large radii and are de¬ 
ficient in gas or show suppressed star formation (e.g., Balogh 
et al. 2000; Solanes et al. 2002; Geha et al. 2012; Wetzel et al. 
2012). Rather, these galaxies may reflect a population of ha¬ 
los that accreted onto their host halo late and are located close 
to the apocenter of their orbit at ~ Rjp- 

As Rsp corresponds to the physical boundary between the 
inner and infall regions of halos, Rsp may also be a natural 
boundary of the one-halo term in halo models of structure for¬ 
mation. Using such a boundary would result in the reclassi¬ 
fication of some isolated halos as subhalos. In Figure 11, we 
explore the magnitude of this effect by comparing the fraction 
of density peaks which are subhalos (i.e., the satellite frac¬ 
tion) as a function of their peak circular velocity, Vpeak^ m the 
standard virial definition to that in the splashback case"^. For 
low values of Vpeak ~ 100 kms“\ the satellite fraction based 
on Rsp shows an increase of 30 percent. 

The reclassification of isolated halos also has implications 
for assembly bias. It is well-known that the large-scale bias of 
halos depends upon their assembly history (Gao et al. 2005; 
Wechsler et al. 2006; Li et al. 2008; Dalai et al. 2008b). At 
the low (high) mass end, the large scale bias of halos is (anti-) 
correlated with their concentrations. The low-mass halos that 
are classified as isolated in the virial definition but as subha¬ 
los in the splashback definition are not a random subsample 
of the parent population. They tend to be in the vicinity of 
more massive halos, have higher concentrations, and higher 
bias values. Thus, their reclassification as satellites can re¬ 
duce the assembly bias of low-mass halos to some extent, but 
it does not explain the bulk of the effect (see also Wang et al. 
2009; Ludlow et al. 2009). 

In halo models, accurate modeling of the transition region 
between the one- and two-halo terms has remained a chal¬ 
lenge. The use of Rsp as the boundary in such models may 
lead to some progress, but requires an accurate calibration of 
the halo mass function and bias as a function of Mjp using 
cosmological simulations. We could convert existing mass 
function and bias formulae for M 2 oom (e.g.. Tinker et al. 2008, 
2010) using the conversions in Equations 7 and 8, but given 
the large scatter in those relations the results would be approx¬ 
imate. Instead, we will require robust methods to measure Rsp 
and Msp for individual halos in simulations, rather than using 
the mean or median density profiles of halo samples. 

4.5. Possible Observational Detections of the Splashback 

Radius 

Interestingly, the steepening of the density profile associ¬ 
ated with the splashback radius may already have been de¬ 
tected in several galaxy clusters. Rines et al. (2013) used the 
velocity caustic method to derive mass profiles of a number of 

^ We have calculated the accretion rate T for individual halos in the Bolshoi 
simulation and assumed zero scatter between T - Ssp in order to define the 
boundaries of these halos. Subhalo classifications were then revised based on 
their updated boundaries. 
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clusters to large radii. Figure 11 in their paper shows that the 
radial density profiles of massive clusters sharply transition to 
slopes steeper than the value of -3 expected from the NFW 
profile. The steepening occurs at radii where p(r) ~ 10-20pm, 
roughly where the splashback radius is expected to be located 
for typical halos with F ~ 1 - 2 (see Figure 2 above and Figure 
2 in Diemer & Kravtsov 2014). 

Furthermore, a deep Chandra observation of Abell 133 
(Vikhlinin et al., in preparation) shows a strong steepening of 
the gas density profile just beyond /? 200 c- Although the steep¬ 
ening is observed in gas rather than in the total mass profile, 
hydrodynamic simulations show that the gas profile exhibits a 
steepening at the same radius as the dark matter profile (Lau 
et al. 2014). According to the calibration of Equation 5, the 
proximity of the splashback radius to /? 200 c would correspond 
to a very high mass accretion rate onto this relatively relaxed, 
“cool core” cluster. In the future, the density profile steepen¬ 
ing associated with splashback may be independently verified 
with mass profile measurements using weak lensing or via ra¬ 
dial profiles of number density of galaxies. 

Finally, Tully (2015) recently argued that a radius qual¬ 
itatively similar to the splashback radius which the author 
called “the second turnaround radius” should be used to define 
the halo extent in analyses of observational samples. Tully 
(2015) presents evidence for a sharp drop in the number den¬ 
sity and velocity dispersion profiles of galaxies in a num¬ 
ber of systems from cluster-sized to the Local Group. We 
note that for the Milky Way and M31, Tully (2015) derives 
a splashback radius of a; 290 kpc, considerably smaller than 
what we would expect for systems of this mass on average; 
Rsp ^ 1.4 - 1.5f?200m ~ 400 - 600 kpc (see Figure 4). Overall, 
Tully (2015) find a scaling of Rgp ~ 1.33f?200c, significantly 
smaller than our calibration. This can be traced to some of the 
simplifying assumptions made in his calculation of the radius. 

4.6. How to Estimate cind M^p in Practice? 

In the standard spherical overdensity definition, the radius 
enclosing a given overdensity A depends on the amount 
of mass within /?a, but not on the distribution of this mass. 
In contrast, estimating M< 4 r, demands knowledge of the scale 
radius (or concentration; Equation 9), and estimating M^p re¬ 
quires a measurement of the density profile at the outermost 
radii of a halo, or knowledge of its accretion rate (Equation 5). 
In simulations, computing M^ 4 r^ presents no problem, since 
the scale radius is routinely computed by halo finders. The 
splashback radius can be found by finding the radius of steep¬ 
est slope in halo density profiles, or by fitting those profiles 
with the fitting function of Diemer & Kravtsov (2014). It may 
be possible to develop more robust methods not relying on 
spherically averaged profiles. Eor example, a clearer picture 
of the density caustic emerges from the distribution of halo 
particles on the phase-space plane (vr - r, where Vi is the ra¬ 
dial velocity with respect to the halo center; Adhikari et al. 
2014). 

Observationally, measurements of r^, R^^p, or the accretion 
rate are challenging, but M< 4 r^ and M^p can be estimated using 
their average relations with conventionally defined masses. 
Eor M^ 4 r^, all we need is a well-calibrated model for halo con¬ 
centrations to estimate r^ (e.g., Diemer & Kravtsov 2015). Eor 
R^p and M^p, we have quantified R^p and M^p as a function of 
f^ 200 m and M 2 oom and the mass accretion rate, given in Equa¬ 
tions (5) and (6). Eor the case where the mass accretion rate of 
a halo is not known, we have given R^p and M^p as a function 
of peak height (or mass) in Equations (7) and (8). However, 


we note that there is one easily accessible quantity that con¬ 
tains information beyond halo mass: the concentration. Since 
the concentration is intimately related to the mass accretion 
history of a halo (e.g., Ludlow et al. 2013), we expect there 
to be a relation between F and concentration. We find that, 
in principle, this relation depends on redshift and mass in a 
non-trivial manner. However, at z = 0 we find a simple, mass- 
independent relation, 

F^ 3.43-2.74 logio(cvir). (10) 

We note that this relation was only calibrated for the cosmol¬ 
ogy used in this paper, and for our particular definition of F. 
Nevertheless, the relation can be used to estimate the mass 
accretion rate, and thus R^p and Mjp, at z = 0. 

For convenience, we have implemented all the calibrations 
of R^p, Msp, and M< 4 r, given in this paper, as a function 
of mass accretion rate, peak height, and mass, in the pub¬ 
lic python code Colossus. This stand-alone module can be 
downloaded at www.benediktdiemer .com/code. 

5. CONCLUSIONS 

We have closely examined the definition of the boundary 
of CDM halos and the associated mass. We have argued that 
the most natural definition of physical halo boundary is the 
splashback radius corresponding to the apocenters of orbits 
of the most recently accreted matter. As shown by Diemer & 
Kravtsov (2014, see also Adhikari et al. 2014), this splashback 
radius manifests itself as a sharp steepening of the outer den¬ 
sity profiles of halos (illustrated in Eigures 1 and 2). In this 
paper, we calibrate the dependence of the splashback radius 
on the halo mass accretion rate and mass (characterized by the 
peak height v). We present a detailed comparison of the evo¬ 
lution of the splashback radius and mass to the evolution of 
radii and masses in the commonly used spherical overdensity 
definitions. In the Appendix we present a detailed analysis 
of the components that contribute to the growth of halo mass. 
Our main findings and conclusions are summarized below. 

1. The splashback radius depends primarily on the mass 
accretion rate of halos, with some cosmological depen¬ 
dence on the mean matter density of the universe, Om(2) 
(Equations (5) and (6)). Eor halos rapidly accreting 
mass, R^p < R^oom, while for slowly accreting halos 
R^p can be as large as ~ 1.5 - 1.6/?200m, up to a; 27?200c- 
Thus, a halo and its environmental effects may extend 
well beyond the commonly used “virial” radii. 

2. A comparison of the growth rates of M^p and masses in 
the standard definitions M 200 C, Mm and Mioom shows 
that, at low z, Msp increases/aifer than M 200 C, Myir, and 
M 2 oomj indicating a substantial physical mass growth at 
low z. Even the smallest halos continue to accrete new 
matter which may fuel their star formation. 

3. To characterize the mass evolution in the inner re¬ 
gions of halos, we introduce a second new mass defi¬ 
nition, M< 4 r,, the mass enclosed within 4 scale radii of 
a halo. This mass is manifestly not affected by pseudo¬ 
evolution because it is not tied to any reference density. 
In the fast mass growth regime, M< 4 r^ =» Myir ~ Mjp, 
while in the slow mass growth regime M< 4 r^ approaches 
a constant. In the latter regime, the interpretation of the 
evolution of the M* - M< 4 r, relation thus becomes very 
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Figure 12. Fractional contribution of pseudo-evolution to halo mass growth, /bdry ■ Left panel: The fraction of mass growth due to the shift in the boundary 
term, between redshift z and z = 0 for halos of different virial masses at z = 0 (Equation (A5)). Center panel: Same as the left panel, but showing the 
instantaneous fraction at redshift z rather than the overall fraction since z (Equation (A6)). Right panel: Same as the center panel, but as a function of Mva- at 
redshift z rather than as a function of the descendant mass. 


simple: at fixed M< 4 r,, the evolution of the relation rep¬ 
resents the change in the stellar mass of a given halo 
along its main progenitor branch. We derive the star 
formation histories of low-mass galaxies in this manner, 
and show that dMtldt{M^ 4 r^ - const) is in good agree¬ 
ment with SFHs inferred from more complex modeling. 

4. A comparison of the evolution of M< 4 r^ and Mjp shows 
that the mass within the inner regions of many halos 
evolves slowly, even though the halo experiences sub¬ 
stantial mass growth overall. The newly accreted mass 
is deposited with a relatively shallow radial density pro¬ 
file and thus contributes most significantly to the outer 
regions. The co-existence of a relatively slow evolu¬ 
tion of the interior mass distribution with an active ac¬ 
cretion in the outskirts means that systems that appear 
quiescent may actually actively grow, particularly be¬ 
cause much of the mass growth is due to the accretion 
of diffuse mass and small halos rather than due to ma¬ 
jor mergers. Thus, in cluster-sized halos, relaxed “cool- 
core” inner regions can co-exist with active accretion in 
the outer regions. 

5. We showed that the relation between the mass accretion 
rates of Mjp and spherical overdensity definitions, M^, 
is complicated by the pseudo-evolution affecting the 
latter. We illustrated the importance of accounting for 
pseudo-evolution by showing that it can accurately ex¬ 
plain the differences in the central and satellite M,-Mvir 
relations found by Watson & Conroy (2013), if halos 
that are satellites at z = 0 were, on average, accreted 
at z = 0.5. After accounting for pseudo-evolution, the 
data imply an even larger amount of post-accretion star 
formation in satellite halos than the 10% increase of 
stellar mass inferred by Watson & Conroy (2013). 

Finally, we alluded to possible observational detections of 
the splashback radius in galaxy clusters. A better observa¬ 
tional characterization of the halo environment close to the 
splashback radius will be an interesting avenue to pursue in 
the near future. Likewise, it will be interesting to explore the 
implications of a larger halo extent of slowly accreting ha¬ 
los for the statistics of their satellite population, the halo mass 


function, bias, and the halo model of galaxy clustering. In par¬ 
ticular, we plan to explore whether statistics such as the halo 
mass function are more universal when recast as a function of 

44sp. 
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APPENDIX 

A. COMPONENTS OE HALO MASS GROWTH 

Here we consider in detail the factors that contribute to 
the evolution of the standard mass definitions employing the 
spherical radius enclosing a given density contrast (see also 
detailed discussion in Section 2 of Diemer et al. 2013). We 
quantify the contribution of mass growth due to the change of 
halo boundary radius for a large range of halo masses and red- 
shifts, and discuss how this contribution should be interpreted 
in the context of halo mass growth. 

Let us consider the change in mass, Ma, of a spherical over¬ 
density halo between two redshifts Zi and Zf. As noted in 
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Diemer et al. (2013), in general this change will consist of 
two components; 

M(zf) = M(zi) + AMtot(zi,Zf) 

— M{z{) + AAfbdry(^i? ^f) ^4^acc(^i; ^f) ? (-^1) 

where 

nzt nRh(z) c/ntr zl 

AMacc(Zi,Zf) = dz I drA-nr^ — (A2) 

Jz, Jo dz 

is due to the actual change of the density profile within the 
halo boundary, while 

rRAizi) 

AMbdry(zi,Zf) = I dr4nr^p(r,Zc)- (A3) 

'^R&(zf) 

is the mass change that arises solely because the boundary of 
the halo increases due to decreasing reference density. We 
have explicitly retained the redshift dependence of the halo 
radius /?a, so that the integral in Equation (A3) contains the 
density at the position r and redshift Zc where Ra(Zc) = r. We 
can further expand this term. 


nK(Z!) 

AMbdry(Zf)= I drAnr"^ p{r,Zi) 

jR(zi) 

jRi 


dr4nr [pir, Zc) - p{r, Zi)] 


R(zi) 

■R(Z!) 


rR(zi: 

jR(zr) 


drA-jTrp(r,Zc) ■ 


(A4) 


The first integral represents the mass change that would arise 
due to the expanding halo boundary if the density profile at Zi 
were to remain unchanged (static), while the other two terms 
correspond to a changing density profile at intermediate red- 
shifts Zc and would be zero for a static profile. Note that we 
integrate the first term only out to ^(zf) which denotes the 
boundary of the halo at Zf, inferred from the density profile at 
redshift Zi. 

The mass change AM^cc clearly corresponds to the newly 
accreted mass within the radius /?a. The AMbdry component, 
however, in general contains both mass that was accreted dur¬ 
ing the interval Az = Zi -Zf and mass that was accreted before, 
at z > Zi, but located at r > R&(zi). The relative fractions in 
these contributions to the mass growth will depend on the ac¬ 
tual choice of the value of A used to define the halo boundary 
and the evolutionary stage of a halo. 

Let us consider the fraction of mass change due to the evo¬ 
lution of the halo boundary: 

fMry(M[Zf],Zf, Zi) = AMbdry/ AM.ot , (A5) 


This fraction is shown in the left panel of Figure 12 forzf = 0 
as a function of Zi and halo mass. For a halo of mass M(zf) 
identified at redshift Zf, the fraction of the instantaneous mass 
growth in its main progenitor at redshift z due to the change 
in its boundary is given by 


Mnst,prog 

'^bdry 


(M[zf],Zf,z) 


liAMbdry) //c/AMtot) 

dz )J [ dz // 


(A6) 


This fraction is shown in the middle panel of Figure 12. 

Finally, the instantaneous fraction of mass growth due to 
the boundary change for a halo of mass M at redshift Zf is 


obtained by taking the limit z ^ Zf of the above equation, 

Zf) = Inn /“:;’'’™*^(M[zf ], Zf , z) (A7) 

and is shown in the right panel of Figure 12. The left and 
central panels can be used to estimate the contribution of the 
boundary evolution to the growth of halo mass along the main 
progenitor branch of a z = 0 halo. The fractions displayed 
in Figure 12 are computed using the mass accretion history 
and concentration model of Zhao et al. (2009). We compute 
the density profile of the main progenitor, p(r, z), assuming 
the Navarro-Frenk-White density profile (Navarro et al. 1997, 
hereafter NFW). 

The left panel of Figure 12 shows that between z = 3 and 
z = 0 the mass change due to the boundary evolution accounts 
for more than half of the change in halo mass, even for mas¬ 
sive 10'"^ /z^'Mq halos. This fraction tends to unity for low- 
mass halos at lower redshifts. Note, in particular, the large 
values of for galaxy-sized halos at z < 1. Moreover, 

> 85% of mass growth for halos of Mvb < 10® Mg, expected 
to host dwarf galaxies, at z < 2 - 3 is due to the halo boundary 
change. 

However, as we discussed in Section 2, the mass increase 
within the splashback radius between two epochs Zi and Zf 
using Equations (A1)-(A3) is due to the new mass shells that 
have entered the splashback radius in the interval Az = Zi - Zf. 
However, AMbdry in Equation (A3) is positive and actually 
represents a significant fraction of the total mass increase, 
AMtot. This illustrates that in the presence of actual mass ac¬ 
cretion, /bdry = AMbdry/AMtot > 0 should not be of particular 
concern. In fact, for the choice of R^^ as the halo boundary, 
the entire AMbdry should correspond to the accretion of new 
mass during the interval Az, and the change of mass within 
should not have any contribution from pseudo-evolution. 
For halo mass definitions using boundaries R < R^^, a frac¬ 
tion of /bdry will correspond to pseudo-evolution rather than 
to new accretion of mass with contribution of PE increasing 
with decreasing physical mass accretion rate. 

In summary. Figure 12 demonstrates that the mass growth 
due to the evolution of the halo boundary is important for all 
halo masses of astrophysical interest, and over a wide range 
of redshifts. How this change should be interpreted, however, 
depends on the evolutionary stage of a halo and on the specific 
choice of A. 
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